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A new numerical algorithm for the double-orbital Hubbard model 
— Hund-coupled pairing symmetry in the doped case 
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In order to numerically study electron correlation effects in multi-orbital systems, we propose 
a new type of discrete transformation for the exchange (Hund's coupling) and pair-hopping in- 
teractions to be used in the dynamical mean field theory + quantum Monte Carlo method. The 
transformation, which is real and exact, turns out to suppress the sign problem in a wide parameter 
region including non-half- filled bands. This enables us to obtain the dominant pairing symmetry in 
the double-orbital Hubbard model, which shows that the spin-triplet, orbital-antisymmetric pairing 
that exploits Hund's coupling is stable in a wide region of the band filling. 

PACS numbers: 71.10.Fd,74.20.Rp,71.30.+h 



In the physics of correlated electron systems, as 
highlighted by the high-T c superconductivity, the pro- 
totypical transition-metal oxides have arrested atten- 
tion because of a rich variety of physical properties, 
among which are anisotropic pairing symmetries in the 
cuprates and in S^RuOj, the colossal magnetoresis- 
tance in manganites, and a complex phase diagram for 
Lai_ x Sr x Mn03. 1] Since the d— orbitals are relevant, we 
are actually talking about multi-orbital systems, and, 
while some of the phenomena should be generically cap- 
tured within single-band models, understanding of the 
effect of electron correlation in multi-orbital system re- 
mains a fundamental problem. 

Namely, the orbital degrees of freedom should couple to 
other degrees of freedom such as charge, spin, and lattice 
distortion, and we can expect even richer physical proper- 
ties. Indeed, the colossal magnetoresistance results from 
Hund's coupling and the complex phase diagrams of Mn 
and Co compounds are considered to arise from a compe- 
tition between Hund's coupling and the Jahn- Teller dis- 
tortion. In perovskite-type oxides, the crystal field splits 
the five 3d— orbitals into three-fold degenerate fag) and 
two-fold degenerate (e g ) levels. When the degenerate 
levels are not fully filled, the degeneracy may be lifted 
by the Jahn- Teller effect, where the system is effectively 
mapped to a single-orbital model when the splitting is 
large enough. We take here the small splitting limit to 
concentrate on the physics specifically caused by the or- 
bital degrees of freedom. 

Both numerical and analytic methods have been de- 
veloped to study correlated electron systems. The meth- 
ods should preferably be non-perturbative if one wants 
to examine the effect specific to the electron correlation 
such as Mott's metal-insulator transition. The dynami- 
cal mean field theory (DMFT)[3j], which can fully include 
temporal fluctuations while spatial fluctuations are ne- 
glected, first succeeded in describing this transition both 
from metallic and insulating sides. In this method a lat- 
tice system (such as the Hubbard model) is mapped to 
an impurity model, which becomes an exact mapping in 



the limit of infinite spatial dimension. A standard proce- 
dure is to solve the impurity problem with the (auxiliary- 
field) quantum Monte Carlo(QMC) method 0, U, which 
involves no approximations except for the Trotter de- 
composition. So the DMFT+QMC method should be 
a desirable candidate for the multi-orbital cases. Un- 
fortunately, it is difficult to extend the QMC method to 
multi-orbitals: First, it is impossible to express those (ex- 
change and pair-hopping) interactions that are specific to 
multi-orbital cases in terms of the usual auxiliary fields. 
Second, even when we can accomplish this, the negative 
sign problem, a notorious problem in QMC calculations, 
is usually difficult to avoid for multi-orbitals. 

This has motivated us to propose here a new auxiliary- 
field transformation that is applicable to Hund's and 
pair-transfer terms. The transformation, which is real 
and discrete, turns out to suppress the sign problem in 
a wide parameter region including non-half-filled bands. 
This enables us to examine the role of Hund's coupling 
and pair-transfer in a double-orbital model. One of 
the most intriguing questions for correlated electrons on 
multi-orbits is what should be the symmetry of the super- 
conducting pairing that arises from the electron-electron 
interaction. We have examined the dominant pairing 
symmetry in the double-orbital Hubbard model with the 
DMFT+QMC method. The result reveals that the spin- 
triplet ® orbital-antisymmetric (g> even-frequency pair- 
ing, which exploits Hund's coupling, is stable in a wide 
region of the band filling. 

So we take the two-fold degenerate Hubbard model[f| 
with a Hamiltonian, 



H — Hq + Hi + Hi , 
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Here cJ m(T creates an electron of spin cr in the orbital 
m(= 1,2) at site i, and n 4mcr = c\ ma c ima . We only con- 
sider the nearest-neighbor hopping between the similar 
orbitals, electron-electron interactions are assumed to be 
intra-atomic with the intra-(inter-)orbital Coulomb in- 
teraction denoted as U(U'), while the exchange and pair- 
hopping interactions as J. The Hamiltonian is rotation- 
ally invariant not only in spin, but also in real space 
if we fulfill the condition U = U' + 2 J (as is the case 
with d-orbitals). We have divided the interaction into 
the density-density interactions Hi and the exchange and 
pair-hopping interactions H 2 . 

The DMFT+QMC method has been used by many 
authors for the single-orbital Hubbard model. For 
the usual on-site Hubbard interaction, a decoupling 
is done with the discrete Hubbard-Stratonovich(HS) 
transformation^ , 



-a[n T rai-|f(n T -fni)] 



lE^ 1 e As(n t-"i)(a > 0) 



±£*V s ("T+'H-i)+f( a <o) 



(2) 



with A = log(e~ + V el°l — 1), which transforms the 
two-body interaction into one-body interactions summed 
over an auxiliary field s. Applying this to each inter- 
action term on the discretized imaginary time, we can 
decompose the partition function Z of the many-body 
system into a sum of the partition functions Z^ s .y of one- 
body systems as Z — J2{ Si } Z{ Si }- The QMC samples 
the single-particle systems according to the weight Z^ s .y 
The negative sign refers to the fact that the weight is not 
positive-definite. The sign problem does not occur, as far 
as the DMFT is concerned, in the single-orbital Hubbard 
model because the impurity problem lacks the electron 
hopping terms. 

While the auxiliary-field QMC method has been ap- 
plied to some multi-orbital Hubbard models by neglect- 
ing the terms other than the density-density interactions 
(Hi in eq.JU), the QMC algorithm becomes a challenging 
problem for the exchange and pair-hopping interactions 
(H2): the HS transformation (j2J is obviously inapplicable 
to these terms. While a decoupling, exp( J Atc\c2c\c4) — 
(1/2) exp[sV JAt(c\c2 — C3C4)], is possible after 
breaking e - ArH ^ into a product of exponentials, it leads 
to a serious sign problem ||. Another attempt by Mo- 
tome and Imada 9] decouples H2 with imaginary aux- 
iliary fields to implement a QMC, but an electron-hole 
symmetry has to be assumed to avoid phase cancellation 
of the weights which become complex due to the imagi- 
nary auxiliary fields. Since the assumption dictates the 
half-filled band and U = U' (that violates U = U' + 2 J), 
investigation of non-half- filled bands and/or the rotation- 
ally symmetric case is difficult with this method. 
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FIG. 1: The average sign plotted as a function of band filling 
at 13= l/T =2 for U = U'= 1, J = 0.5 and at /3 = 10 for U = 
U =1, J = 0.2 calculated with the present algorithm (circles) 
and with the one due to :9] (squares). In this particular plot 
we have set U = U' to facilitate comparison with Motome and 
Imada's method. [Tl| Curves are guides to the eye here and in 
subsequent figures. 



We propose here a new type of discrete transformation 
for the exchange and pair-hopping interaction terms, 



-AtH 2 



- V e Ar(/ T- / i ) e a(Ar T+ Ar i)+ &Ar T^i 
2 

r=±l 



(3) 



where 



A = ^\og(e 2JAT + Ve iJAT -1), 
a = — log(cosh(A)), b = log(cosh( J At)), 
f a = c\ a C2a + 4 CT c i<r> N a = n la + n 2(T - 2n 1(T ri2<r. 

This transformation is exact with the auxiliary field (r) 
being real and all the operators being hermitian. Al- 
though a term N^N^ which is forth order in n appears 
on the right hand side, we can apply the usual HS trans- 
formation to this term due to a property N% = N a . The 
resulting terms in the form nn can be combined with the 
Coulomb terms, which are transformed with eq.©, so 
that we need in total only two auxiliary-fields for H2. 
Furthermore the interaction parameters U,U' ', J can be 
varied independently, which means that we can treat the 
rotational symmetric cases of U = U' + 2J. We note that 
recently Han|ld| proposed another type of HS transfor- 
mation for the Hamiltonian Q . He adopts a continuous 
HS transformation in combination with the discrete one 
l|2|l. and avoided the sign problem in a wide parameter 
region at half-filling. In contrast all the transformation 
are discrete in the present method. 

One valuable property of the transformation @ is 
that negative weights are indeed reduced significantly. 
An unexpected asset is that the suppression of the sign 
problem is almost independent of the band filling. This 
is shown in Fig^ where the average sign (sign) = 
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FIG. 2: The average sign plotted as a function of inverse 
temperature for n = 1.8, U = U' = 1, J = 0.2 calculated with 
the present algorithm (circles) and with the one due to 0] 
(squares). 



2{ Si } ■^{s i }/X){s i } l^{si}l) 1S plotted as a function of the 
band filling for two sets of values of j3= 1/T and J for the 
two-orbital Hubbard model on an infinite dimensional hy- 
percubic lattice, whose density of states is Gaussian with 
the band width W = 2. We can see that the sign problem 
is avoided almost independently of the band filling. Fig- 
ure [3 depicts the temperature dependence, which shows 
that (sign) decreases mildly as the temperature lowers, so 
we can go to lower temperatures with the present trans- 
formation. One reason why the transformation J3J) re- 
duces the amount of negative weights is that it has a 
single auxiliary field (r in eq.©) as the source of the 
negative sign; the other auxiliary fields, related to the 
density-density interactions, are irrelevant to the nega- 
tive weights. 

To test how the present method works, we have applied 
it to the calculation of the superconducting susceptibility 
for the Hamiltonian Since off-site, anisotropic pair- 
ings (such as p and d- waves) cannot be treated within the 
DMFT, we confined ourselves to s-wave pairing. Even 
within that channel, various pairing symmetries are pos- 
sible in multi-orbital systems, since the total symmetry 
now consists of spin (g> orbital ® frequency, so 
spin orbital frequency 



ISE 
3AE 
IAO 
3SO 



singlet symmetric even 

triplet antisymmetric even 

singlet antisymmetric odd 

triplet symmetric odd 



are the possibilities. The pairs that are formed across 
different orbitals are especially interesting. The orbital- 
symmetric pairs are 



S a : ci T c u + c 2T c 2i , 
S b : ci T c u - c 2T c 2 j., 
S c : ci T c 2 j, + C2TC1J., 



(4) 



where cifCij and C2fC 2 x are combined into bonding and 
antibonding states due to the pair-hopping term, while 
the orbital-antisymmetric pair is 



A 



Ci T C 2 | - C2TC1J.. 



Here we take, without the loss of generality, the S z = 
channel for spin triplets, where S z is the z-component of 
the Cooperon spin. 

So the interest here is which pairing symmetry is fa- 
vored in the double-orbital system, especially in the pres- 
ence of Hund's coupling. For the single-orbital case the 
sites are mostly singly-occupied by electrons around the 
half-filling, where the pairing (usually in d-wave channel) 
occurs. The question is what would be the corresponding 
picture for the double-orbital case around half-filling. So 
we have calculated the pairing susceptibility P, which is 
related to the two-particle Green function|12|. 

(iw n ,iaj' n ) 
pP r p r p r p 

Jo Jo Jo Jo 

x E( T ^fc/' T ( r i) c -fc/i( T 2) ct _ fc ' m i( T3 ) c fc' m 'T (r4)> (5) 

kk' 



through the equation, 



P=g^xg= Yl 9u>{u n )xil> 

W inm' uj n u)' 



,(iu) n ,iu' n )g mm ,{uj' n ). 



Here g is the form factor describing the pairing symmetry 
in orbital-frequency space, which is either even or odd 
in w for each orbital component. We adopt <7odd(^) = 
sign(w) and g C vcn(^) = 1 here. 

The two-particle Green function © is obtained from 
the Bethe-Salpeter equation, 



X = Xo + Xo^x, 



(6) 



where \ is defined in eq.(|SJ, Xo the bare two-particle 
Green function calculated as a product of the one-particle 
Green function, G|(fc, iw)Gj,(— fc, — iui), summed over the 
momenta, and T the vertex function, all of them being 
matrices with respect to orbital and frequency indices. 
Within the DMFT we can replace T with the local one, 



^loc 



= (xrTr 1 - (x 1o T\ 



in the limit of infinite dimensions ^3], where Xo° c an d 
X loc are respectively the bare and the dressed two-particle 
Green functions for the effective impurity model, which 
are computed in the QMC. 

The temperature dependence of the pairing suscepti- 
bilities is shown in FigCflfor n = 1.8, V = 0.7, J = 0.4, 
and U = U'+2J= 1.5. We can see that the spin-triplet, 
orbital-antisymmetric, even frequency pairing (denoted 
as 3AE) becomes dominant at low temperatures. The 
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FIG. 3: Temperature dependence of the pairing susceptibil- 
ities for n = 1.8, U' = 0.7, J = 0.4, U = U' + 2J = 1.5. 
The symbols denote pairing symmetries (1: spin-singlet, 3: 
triplet; S a,b,c : orbital-symmetric (eq.@), A: antisymmetric; 
E: even-, O: odd- frequency) . 
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FIG. 4: Band-filling dependence of the pairing susceptibilities 
for /3 = 60, U'=0.7, J = 0.4, U = U' + 2 J = 1.5 with the same 
abbreviation for pairing symmetries as in Fig[3] 



enhancement of such a pairing should be due to Hund's 
coupling, which tends to align electron spins across dif- 
ferent orbitals. The result is consistent with Han's[l"(1] 
for the semi-elliptical density of states. While in the 
work of Han the susceptibility diverges at T/W ~ 0.06 
for J/W = 0.15, U'/W=0A5, we do not detect the diver- 
gence up to T/W = 1/180. This should be because elec- 
trons are less correlated in the hypercubic lattice than in 
the Bethe lattice (which the semi-elliptic density of states 
would represent), since the Gaussian density of states for 
the hypercubic lattice, with high-energy tails, has an ef- 
fectively larger band width. 

Figure|Hsliows the band-filling dependence of the pair- 
ing susceptibilities at T/W= 1/120 for 17' = 0.7, J = 0.4, 



and U = U'+2 J=1.5. A new finding here is that the spin- 
triplet, orbital-antisymmetric (3AE) pairing induced by 
Hund's coupling remains dominant in a rather wide range 
(n = 1.6 ~ 2) of hole-doping at this temperature. It de- 
clines, however, at low fillings (n = 1 ~ 1.4). This re- 
sult may be related to a mechanism of superconductivity 
proposed by Capone0| for multi-orbital systems close 
to the Mott transition, where the electron repulsion U is 
envisaged to assist a pairing even for s- waves. 

In summary, we have constructed a new discrete trans- 
formation (which reduces the sign problem and is appli- 
cable to doped bands) for the exchange and pair-hopping 
terms in double-orbit models, and implemented this in a 
DMFT+QMC calculation. Superconducting susceptibil- 
ities of the s-wave pairings calculated with this method 
shows that the spin-triplet, orbital- antisymmetric even 
frequency paring, enhanced by Hund's coupling and be- 
ing dominant at half-filling, is robust against hole doping. 
We note that it is difficult to investigate ferromagnetism 
in this model even with the present transformation, since 
the ferroma gnet ism appears only in very strongly corre- 
lated regime [lj|, where the sign problem becomes serious 
again, so this is a future problem. 

Numerical calculations were partly performed on 
SR8000 in ISSP, University of Tokyo. 
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